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I Abstract 

Clear signs of excess dileptons above the known sources were found at the SPS since long. 

I However, a real clarification of these observations was only recently achieved by NA60, measuring 

dimuons with unprecedented precision in 158A GeV In-In collisions. The excess mass spectrum 
in the region M<1 GeV is consistent with a dominant contribution from tt+tt" ^ p — > ii~ 
annihilation. The associated p spectral function shows a strong broadening, but essentially no shift 

I in mass. In the region M>1 GeV, the excess is found to be prompt, not due to enhanced charm 

production. The inverse slope parameter Tefj associated with the transverse momentum spectra 
rises with mass up to the p, followed by a sudden decline above. While the initial rise, coupled to a 

\ hierarchy in hadron freeze-out, points to radial flow of a hadronic decay source, the decline above 

signals a transition to a low-flow source, presumably of partonic origin. The mass spectra show 
the steep rise towards low masses characteristic for Planck-like radiation. The polarization of the 

\ excess referred to the Collins Soper frame is found to be isotropic. All observations are consistent 

with a global interpretation of the excess as thermal radiation. We conclude with a short discussion 
of a possible link to direct photons. 

\ Introduction 

Dileptons are particularly attractive to study the hot and dense QCD matter formed in high-energy 
. nuclear collisions. In contrast to hadrons, they directly probe the entire space-time evolution of the 
' expanding system, escaping freely without final-state interactions. At low masses M<1 GeV (LMR), 
\ thermal dilepton production is mediated by the broad vector meson p (770) in the hadronic phase. 
■ Due to its strong coupling to the vrvr channel and the short life time of only 1.3 fm/c, "in-medium" 
modifications of its mass and width close to the QCD phase boundary have since long been considered 
as the prime signature for chiral symmetry restoration [H [21 [3]. At intermediate masses M>1 GeV 
(IMR), it has been controversial up to today whether thermal dileptons are dominantly produced in the 
earlier partonic or in the hadronic phase, based here on hadronic processes other than vrvr annihilation. 
Originally, thermal emission from the early phase was considered as a prime probe of deconfinement [HE]. 

Excess dileptons above the known decay sources at SPS energies were observed before by CERES 
[61 Ej for M<1 GeV, NA38/NA50 [8] for M>1 GeV and by HELIOS-3 [9J for both mass regions 
(see [To] for a short recent review including the preceding pp era and the theoretical milestones). The 
sole existence of an excess gave a strong boost to theory, with hundreds of publications. In the LMR 
region, vrvr annihilation with regeneration and strong in-medium modifications of the intermediate p 
during the fireball expansion emerged as the dominant source. However, the data quality in terms of 
statistics and mass resolution remained largely insufficient for a precise assessment for the in-medium 
spectral properties of the p. In the IMR region, thermal sources or enhanced charm production could 
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Figure 1: Left: Background-subtracted mass spectrum before (dots) and after subtraction of the known 
decay sources (triangles). Right: Excess dimuons compared to theoretical predictions [17], renormalized 
to the data in the mass interval M<0.9 GeV. No acceptance correction applied. 



account for the excess equally well, but that ambiguity could not be resolved, nor could the nature of 
the thermal sources be clarified. 

A big step forward in technology, leading to completely new standards of the data quality in this field, 
has recently been achieved by NA60, a third-generation experiment built specifically to follow up the 
open issues addressed above [T^. Results on mass and transverse momentum spectra of excess dimuons 
have already been published [TT|, [121 [IS]- Further aspects associated with the centrality dependence, 
polarization and acceptance-corrected mass spectra were reported during 2008 [HI [15]. This paper 
shortly reviews the central results and also discusses a possible link of the excess dileptons at very low 
masses to real photons. 

Mass and spectra of excess dileptons 

Fig. [1] (left) shows the centrality-integrated net dimuon mass spectrum for 158A GeV In- In collisions 
in the LMR region. The narrow vector mesons u and are completely resolved; the mass resolution at 
the u is 20 MeV. The peripheral data can fully be described by the electromagnetic decays of neutral 
mesons p!2| [16] . This is not true for the more central data as plotted in Fig. [H due to the existence 
of a strong excess of pairs. The high data quality of NA60 allows to isolate this excess with a priori 
unknown characteristics without any fits: the cocktail of the decay sources is subtracted from the total 
data using local criteria, which are solely based on the mass distribution itself. The p is not subtracted. 
The excess resulting from this difference formation is illustrated in the same figure (see [121 131 [IS] for 
details and error discussion). 

The common features of the excess mass spectra can be recognized in Fig. [1] (right). A peaked 
structure is always seen, residing on a broad continuum with a yield strongly increasing with centrality 
(see below), but remaining essentially centered around the nominal p pole [TB]. Without any acceptance 
correction and px selection, the data can directly be interpreted as the space-time averaged spectral 
function of the p, due to a fortuitous cancellation of the mass and pt dependence of the acceptance 
filtering and the phase space factors associated with thermal dilepton emission [16]. The two main 
theoretical scenarios for the in- medium spectral properties of the p, broadening [2] and dropping mass [3], 
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are shown for comparison [T7j. Clearly, the broadening scenario gets close, while the dropping mass 
scenario in the version which described the CERES data reasonably well [21 El E] fails for the much 
more precise NA60 data. 
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Figure 2: Excess yield ratios for peak, con- 
tinuum and total vs. centrality for the mass 
window 0.2<M<1 GeV and all p^. 
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A detailed view of the shape of the excess mass spectra is obtained by using a side window 
method [16] to determine separately the yields of the peak and the underlying continuum. Fig. [2] 
shows the centrality dependence of these variables: peak, underlying continuum and total yield in the 
mass interval 0.2<M<1.0 GeV, all normalized to the cocktail p. The continuum and the total show 
a very strong increase, starting already in the peripheral region, while the peak slowly decreases from 
>1 to <1. Recalling that Fig. [T] roughly represents the full p spectral function, the excess/p ratio can 
directly be interpreted as the number of p generations created by formation and decay during the fireball 
evolution, including freeze-out: the 'p-clock', frequently discussed in the past. It reaches up to about 6 
generations for central In-In collisions; selecting lower px this number doubles. 
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Figure 3: Left: Fit of the weighted offset distribution in the IMR region with the contributions from 
charm and prompt decays. Right: Acceptance-corrected mass spectra of Drell-Yan, open charm and 
the excess (triangles). 



The central NA60 results in the IMR region [TT] are shown in Fig. [31 The use of the S'i-vertex 
tracker allows to measure the offset between the muon tracks and the main interaction vertex and 
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Figure 4: Acceptance-corrected transverse mass spectra of the excess dimuons for 4 mass windows and 
the (p (left), and a decomposition into peak and continuum for the p-hke window (right, see text). 
The hues represent exponential fits to the data in the range 0.4<pt<1.8 GeV. The normalization in 
absolute terms is independent of rapidity over the region measured. For error discussion see [13]. 



thereby to disentangle prompt and offset dimuons from D decays. The left panel of Fig. [3] shows 
the offset distribution to be perfectly consistent with no charm enhancement, expressed by a fraction 
of 1.1 ±0.1 of the canonical level. The observed excess is really prompt, with an enhancement over 
Drell-Yan by a factor of 2.4. The excess can now be isolated in the same way as was done in the 
LMR region, subtracting the measured known sources, here DY and open charm, from the total data. 
The right panel of Fig. [3] shows the decomposition of the total into DY, open charm and the prompt 
excess. The mass spectrum of the excess is quite similar to the shape of open charm and much steeper 
than DY; this explains of course why NA50 could describe the excess as enhanced open charm. The 
transverse momentum spectra are also much steeper than DY [H]. The fit temperatures of the ttlt 
spectra associated with 3 mass windows are indicated on the bottom of the figure. 

The remainder of this paper is solely concerned with excess data fully corrected for acceptance and 
pair efficiencies (T3t [18] . In principle, the correction requires a 4-dimensional grid in the space of M- 
PT-y-cos9cs (where 9cs is the polar angle of the muons in the Collins Soper frame). To avoid large 
statistical errors in low-acceptance bins, it is performed instead in 2-dimensional M-pj- space, using 
the measured y and cos6 distributions as an input. The latter are, in turn, obtained with acceptance 
corrections determined in an iterative way from MC simulations matched to the data in M and pt- 
The ^/-distribution is found to have the same rapidity width as dNch/dr], ay ~ 1.5 [18]. The cosOcs 
distributions for two mass windows of the excess and the uj are contained in |T3]. Within errors, they 
are found to be uniform, implying the polarization of the excess dimuons to be zero, in contrast to DY 
and consistent with the expectations for thermal radiation from a randomized system. 

The two major variables characterizing dileptons are M and pt-, and the existence of two rather than 
one variable as in case of real photons leads to much richer information. Beyond (minor) contributions 
from the spectral function, pt encodes the key properties of the expanding fireball, temperature and 
transverse (radial) flow. In contrast to hadrons, however, which receive the full asymptotic flow at the 
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Figure 5: Left: Acceptance-corrected mass spectra of the excess dimuons for the combined LMR/IMR 
regions. The theoretical scenarios are labeled according to the authors RR [19], HR [20], and DZ |21j . 
In case of [20], the EoS-B"*" option is used. Right: Inverse slope parameter Tgfj vs. dimuon mass for the 
combined LMR/IMR regions of the excess in comparison to hadrons [13]. Open charm is subtracted in 
both parts of the figure. 



moment of decoupling, dileptons are continuously emitted during the evolution, sensing the space-time 
development of temperature and flow. This makes the dilepton px spectra sensitive to the emission 
region, providing a powerful diagnostic tool [5l[l9]. Fig. H] (left) displays the centrality-integrated uit 
spectra, where my = (p^+M2)V2^ for four mass windows; the is included for comparison. The ordinate 
is normalized to dNch/drj in absolute terms [HJ [15]. Apart from a peculiar rise at low rriT (<0.2 GeV) 
for the excess spectra (not the 0) which only disappears for very peripheral collisions [TOl [T3] , all spectra 
are pure exponentials, but with a mass- dependent slope. Fig. |4] (right) shows a more detailed view into 
the p-like mass window, using the same side-window method as described in connection with Fig. [2] 
to determine the pt spectra separately for the p peak and the underlying continuum. All spectra are 
purely exponential up to the cut-off at pr=3 GeV, without any signs of an upward bend characteristic 
for the onset of hard processes. Their slopes are, however, quite different (see below). 

The inverse slope parameters Tes extracted from exponential fits to the ttlt spectra (with a finer 
binning than in Fig. Hj) are plotted in Fig. [5] (right) vs. dimuon mass, unifying the data from the LMR 
and IMR regions. The hadron data for rj, uj and obtained as a by-product of the cocktail subtraction 
are also included, as is the single value for the p-peak from Fig. H] (right). Interpreting the latter as 
the freeze-out p without in-medium effects, all four hadron values together with preliminary vr" data 
from NA60 can be subjected to a simple blast- wave analysis [2]. This results in a reasonable set of 
freeze-out parameters of the fireball evolution and suggests the following consistent interpretation of 
the hadron and dimuon data together. Maximal flow is reached by the p, due to its maximal coupling 
to pions, while all other hadrons freeze out earlier. The Teg values of the dilepton excess rise nearly 
linearly with mass up to the p-pole position, but stay always well below the p line, exactly what would 
be expected for radial flow of an in-medium hadron-like source (here Tr+vr^^p) decaying continuously 
into dileptons. 

Beyond the pole region of the p, however, the Tes values of the excess show a sudden dechne by about 
50 MeV. Extrapolating the lower-mass trend to beyond the p, such a fast transition to a seeming low- 
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flow situation is extremely hard to reconcile with emission sources which continue to be of dominantly 
hadronic origin in this region. A more natural explanation would then be a transition to a dominantly 
early partonic source with processes like qq Z^"*"/^" for which flow has not yet built up [12]. While 
still controversial pO] , this may well represent the first direct evidence for thermal radiation of partonic 
origin, overcoming parton-hadron duality for the yield description in the mass domain. 

The acceptance- and efficiency-corrected data can also be projected on the mass axis. The p^- 
differential results are discussed in [HI [15]. Fig. [5] (left) shows pT-integrated data for pt>0.2 GeV (the 
lowest pt bin has been excluded to avoid the very large errors in the low-acceptance region [TU [15]). 
The LMR and IMR regions are again unified, using an independently assessed absolute normaliza- 
tion [22]. Recent theoretical results from the three major groups working in the field are included for 
comparison [T9|, [20|, [21] . A strong rise towards low masses is seen in the data, reflecting the Boltzmann 
factor, i.e. the Plank-like radiation associated with a very broad, nearly fiat spectral function. Only 
the Hees/Rapp scenario [2D] is able to describe this part quantitatively, due to their particularly large 
contribution from baryonic interactions to the low-mass tail of the p spectral function. The p pole 
remains visible, due to the freeze-out part. In-medium broadening of parts of the u and are contained 
in [20], but not in [l9l[2T], accounting for the broad bump seen in the region of the 0. Caution should, 
however, be taken as to the nearly quantitative agreement with the data, since some details of the 
resolution function of the NA60 apparatus are still under investigation. 



Real photons and low-mass dileptons 

Processes contributing to dilepton production at very low masses include vrvr annihilation (extending 
to below 2Mt^ because of the softening of the pion dispersion relation), vr-baryon processes, the ai-Dalitz 
decay, parton annihilation, internal conversion of photons from processes creating also real photons 
(e.g. Compton-like graphs or bremsstrahlung) , and others [Hj, [201 [21], and the superposition of all 
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Figure 6: Acceptance-corrected pr 
spectrum of the excess for the lowest 
mass bin 0.2<M<0.3 GeV. The solid 
line corresponds to an exponential fit 
(as in Fig. [1]) for pt>1 GeV, while 
the dashed line is just an extrapolation 
down to Pt<1 GeV. 



these sources accounts for the total measured yield as shown in Figs. [T] and [5] PHENIX at RHIC has 
recently [23] analysed dielectron data in the mass window 0.2<M<0.3 GeV for Pt>1 GeV under the 
extreme assumption, that the total yield for that selection would solely be due to internal conversion 
of real photons; the transformation back to real photons was done with the prescription in [2l] . The 
analogous situation for NA60 is illustrated in Fig. [6], but with the complete pr-spectrum; it corresponds 
to the lowest point in Fig. [5] (right). The fit value for the cut Pt>1 GeV is T=183±10 MeV. However, 
this value of T corresponds to the effective temperature Teg as used in Fig. [5l affected by radial flow 
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(ignored by PHENIX). If corrected for that, using the intercept of the hnear rise with the ordinate, it 
is reduced to about 150 MeV. Such a low value hardly makes sense in the light of initial temperatures 
of 210-350 MeV expected for real photons [25], but is more indicative for the average temperature of 
the hadronic phase alone. In Fig. [5l the lowest point does not look exceptional at all, but just appears 
as part of the general systematic rise of Tes vs. mass. The lack of any convincing tag on internal 
conversion and the embeddment of the lowest mass bin in the general systematics rule out the real 
photon hypothesis for this bin, and this will hardly change at higher masses. 

Conclusions 

Excess dileptons have been observed over the complete mass region 0.2<M<2.6 GeV. Their p^ 
spectra offer a clear view of the origin of the different sources. All observations, including the finding of 
zero polarization, are consistent with a global interpretation of the excess as thermal radiation. There 
is no way to unambiguously isolate real photons via internal conversion of low-mass dileptons. 
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